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ABSTRACT 

"Afternoon  effect"  is  a  type  of  sonar  performance  loss  encountered 
during  the  afternoons  and  evenings  of  calm  sunny  days.  It  is  caused  by 
a  temperature  rise  near  the  top  of  the  otherwise  isothermal  mixed  layer, 
increasing  the  sound  speed  there  and  so  refracting  the  sound  rays 
downwards,  away  from  the  sonar  target.  The  temperature  rise  profile  is 
governed  by  wind  speed,  solar  heating  and  other  meteorological  factors, 
so  that  afternoon  effect  prediction  should  become  an  extension  of  meteor¬ 
ological  forecasting.  However  there  is  at  present  no  established  method 
for  computing  the  temperature  profile  from  the  meteorological  information. 

A  mixed-layer  model  is  presented  which  appears  to  meet  this  need, 
subject  to  further  testing.  The  influences  of  some  key  environmental 
factors  are  explored.  An  example  of  afternoon  effect  prediction  is  given, 
and  the  next  steps  towards  a  routine  prediction  facility  are  indicated. 
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SUMMARY 

"Afternoon  effect"  is  a  type  of  sonar  performance  loss 
encountered  during  the  afternoons  and  early  evenings-  of  calm  sunnv  days. 

It  is  caused  by  the  development  of  a  temperature  rise  near  the  top  of  the 
otherwise  mixed  layer,  known  as  a  "transient  thermot line"  or  simply 
"transient".  This  increases  the  local  sound  speed,  thus  refracting  sound 
waves  downwards  away  from  the  target.  In  round  terms  the  loss  is  or  is  not 
significant,  depending  on  whether  the  sound  speed  at  sonar  projector  depth 
does  or  does  not  exceed  the  speed  at  the  bottom  of  the  mixed  laver.  'Hie 
mixed  layer  sound  speed  profile  depends  upon  the  temperature  profile,  which 
in  turn  is  governed  by  wind  speed,  solar  heating  and  other  meteorologies , 
factors.  Hence  the  prediction  of  afternoon  effect  should  become  an 
extension  of  conventional  meteorological  forecasting. 

However  there  is  as  yet  no  established  method  for  determining  the 
diurnal  temperature  profiles  from  the  meteorological  data.  A  model  is  now- 
presented  which  appears  to  meet  this  need  subject  to  some  further  testing. 

The  model  uses  the  standard  one-dimensional  eddy-dif fusivitv  f K ) 
equations  for  heat  and  horizontal  momentum.  These  are  closed  by  a  third 
dif fusivity-type  equation  in  K  itself,  the  generation  term  of  which  derives 
from  the  standard  balance  for  turbulent  kinetic  energy.  The  last  contains 
a  prescription  for  turbulence  dissipation,  including  a  turbulence  length 
scale,  adapted  from  the  literature.  The  three  diffusion  equations  are 
solved  by  well-known  numerical  methods  with  computer  assistance.  The 
solutions  are  well-behaved  and  do  not  require  smoothing. 

The  model  is  fitted  to  a  data  set  acquired  by  the  author,  and 
then  tested  against  established  empirical  rules  for  transient  formation. 


K  is  found  to  decrease  sharply  with  depth,  through  up  to  five  decades,  to 
a  minimum  at  the  transient  bottom.  The  effect  is  to  restrict  sub-transient 
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diffusion  of  momentum  and  heat,  totally  at  low  wind  speeds;  the  effect 
diminishes  as  the  speed  increases.  This  explains  the  "half-section  of 
a  wineglass"  appearance  characteristic  of  developed  transients.  The 
dissipation  term  in  the  turbulent  energy  balance  is  found  to  be 
comparatively  minor,  connoting  severe  stratification  within  the  transient 
with  gradient  Richardson's  number  Ri  exceeding  one.  This  sets  the  model 
apart  from  others  wherein  Ri  is  limited  to  about  0.25. 

The  influences  of  a  number  of  environmental  factors  are 
explored.  Generally  it  is  concluded  that  any  factor  which  alters  the 
balance  between  the  effects  of  wind  and  solar  radiation  will  have  a 
major  influence  on  transient  formation. 

An  example  of  transient  prediction  is  given,  and  t le  next  steps 


towards  a  routine  afternoon-effect  prediction  facility  are  indicated. 
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FIGURES  (in  text) 


Page  No. 

1.  Some  observed  temperature  profiles  in  the  nominally-mixed  1 

ocean  layer,  at  1600  local  time  except  as  indicated. 

Other  bracketed  figures  are  wind  speed  mean  or  range. 

2.  Profile  obtained  by  special  slow-falling  XBT  from  DTV  SEAL  3 

on  6  Apr  1982  at  1300  local  time,  50nm  S.E.  of  Sydney. 

Wind  speed  2-5  kn,  cloud  cover  1-2  eighths. 

3.  Transients  corresponding  roughly  to  particular  windspeeds,  5 

with  corresponding  sound  velocity  profiles,  and 

hypothetical  sonar  ranges  in  terms  of  receiver  depth. 

4.  Mixed-layer  temperature  differential  for  zero  sound  speed  excess.  7 

5.  HMAS  DIAMANTINA  cruise  May  1976,  temperature  profiles  obtained  by  21 

XBT.  Base  temperatures  have  been  adjusted  to  give  overall  heat 
balance . 

6.  UMAS  DIAMANTINA  cruise  May  4/5  1976,  composited  meteorological  25 

data  :  wind  speed  and  direction,  and  inferred  rates  of  radiant 
heating  and  surface  cooling. 

7.  Radiant  penetration  -  depth  factors  (RZ)  for  representative  29 

oceanic  and  coastal  type  waters. 

8.  HMAS  DIAMANTINA  cruise  May  4/5,  1976.  Upper  panel  :  composited  31 

XBT  temperature  profiles  (heavy  line)  at  1200,  1400,  1600  and 

1800  local  times,  from  Fig,  5;  also  corresponding  modelled 
profiles  (light  line)  with  dissipation  index  P  =  4  and  F- 
values  as  marked,  using  meteorological  data  of  Fig  6.  Water 
assumed  to  be  of  Oceanic  lb  classification.  Lower  panel  : 

The  same  except  that  the  modelled  profiles  are  for 
P  =  2,  and  0,  with  F  optimised  at  values  1600,  4  and 
0.4  respectively. 

9.  Modelled  results  for  a  midsummer  day  in  latitude  53°N,  for  37 

wind  speeds  of  5,  10  and  15  knot.  DNOM  predictions  of  depths 

to  top  and  bottom  of  transient  are  superimposed  in  broken  line 
and  shaded  band  respectively.  The  right-hand  column  shows  the 
eddy-dif f usivity  (K)  results. 

10.  Variation  of  VG2  (square  of  velocity  gradient),  BG  (buoyancy  41 

gradient)  and  DIS  (  =  e/K,  modified  dissipation)  with  depth. 

Conditions  are  as  for  Fig.  9  with  F“  4 

11.  Variations  of  ^  and  £  ff  (K~j  with  depth  for  two  wind  speeds,  ^5 
marked  t  and  z  respectively.  Conditions  are  otherwise  as  for 

Figs.  9  and  10. 

12.  Temperature  profiles  for  local  time  1600  repeated  from  Fig.  9  with  47 
F  ■  4,-  showing  the  effect  of  varying  the  assumed  sunrise 
diffusivity  by  factors  of  0.25  and  4. 

The  diffusion  equations  for  momentum  heat  and  diffusivity  set 
out  in  flowsheet  form,  corresponding  to  the  scheme  of 
numerical  solution  described  in  Section  2e . 


13. 


49 
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14.  Ri-profiles  for  three  wind  speeds  at  rime  1600.  Conditions  5? 

are  as  for  Fig.  9,  with  F  =  4  . 

15.  HMAS  DIAMANTINA  cruise  May  4/5  1976,  temperature  and  Ri  55 

profiles  at  1600  local  time  predicted  by  Mel  lor  and 

Durbin  (1975)  model  formulation.  The  central  a-value 
0*1  is  that  recommended  by  MD.  The  observed  and  nresently- 
modelled  (F  =  4)  profiles  are  repeated  from  Fig.  8  for 
comparison. 

16.  Temperature  profiles  for  a  midsummer  day  in  latitude  3  3N  at  5? 

local  time  1600,  with  wind  speeds  5  and  10  knot,  obtained 

using  Mellor  and  Durbin  (1975)  model  formulation.  The 
central  a-value  0*1  is  that  recommended  bv  MD.  The 
presently-modelled  profile  (  F  =  4 )  and  the  DNOM  predicted 
top  depth  are  repeated  from  Fig.  9  for  comparison. 

17.  Time  1600  temperature  and  diffusivitv  results  repeated  from  59 

Fig.  9  (curves  S) .  The  other  curves  are  for  identical 
conditions  except  for  environmental  changes  :  CC  greatly 
increased  cloud  cover;  TW  turbid  water;  WW  warm  tropical 

water;  NC  no  Coriolis  acceleration. 

18.  "Sydney  exercise"  computed  results.  Conditions  :  Lat .  34°S,  61 

bright  day,  local  time  1600. 
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1.  HMAS  DIAMANTINA  cruise  1-5  .May  1976.  Principal  daytime  features.  24 

2.  Unabsorbed  percentage  of  radiant  flux  at  particular  depths,  for  28 
representative  oceanic  and  coastal  type  waters. 

3.  Transient-depth/wind-speed  relationship  after  Tabata  et  al  (1965,  35 

Fig.  1) 
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Fig.  1.  Some  observed  temperature  profiles  in  the  nominally-mixed 
ocean  layer,  at  1600  local  time  except  as  indicated. 

Other  bracketed  figures  are  wind  speed  mean  or  range. 


1 .  INTRODUCTION 
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1.1  Transient  thermoclines 

"Afternoon  effect"  is  a  loss  of  performance  often  suffered  bv 
surface-duct  sonars  under  conditions  of  bright  sunlight  and  low  wind, 
so  called  because  it  usually  peaks  at  that  time  of  day.  It  is  caused  bv 
the  development  of  a  temperature  gradient  close  to  the  top  of  the  mixed 
layer,  representing  the  accumulation  of  solar  heat  in  the  absence  of 
effective  downward  mixing  by  the  wind.  This  temperature  gradient  produces 
a  sound  velocity  gradient  which  refracts  sound  rays  downward;  in 
consequence  the  latter  can  escape  from  the  surface  duct  altogether,  causing 
the  sonar  performance  loss  already  mentioned. 

These  near-surface  temperature  gradients  are  termed  "transient 
thermoclines",  often  abbreviated  simply  to  "transients",  to  distinguish 
them  from  the  seasonal  thermoeline  lying  beneath  the  mixed  layer.  Transient 
observation  is  well-established  in  the  literature.  Sverdrup  (1942,  Table  32) 
gives  diurnal  surface  temperature  rises  for  tropical  oceans,  typically 
1.6°C  for  a  nil  to  light  breeze  and  0.7°C  for  a  moderate  to  fresh  breeze. 
These  figures  are  for  a  clear  sky;  for  an  overcast  sky  they  fall  to  0 ,9°C 
and  0.4°C  respectively.  More  recently,  sets  of  temperature  profiles  taken 
during  a  single  day  have  appeared  in  the  literature,  see  for  example  Howe 
and  Tait  (1969),  Kaiser  (1978),  Shonting  (1964),  and  Stommel  et  al  (1969). 
Fig.  1  shows  a  late-af ternoon  profile  from  each  of  these  sources,  marked 
with  representative  wind  speed.  A  distinctive  "half-section  of  a  wineglass" 
appearance  is  seen  to  be  common  to  all.  To  these  may  be  added  an  isolated 
profile  acquired  by  the  author  about  50  nro  S.E.  of  Sydney,  in  April  1982  at 
1300  local  time.  This  is  shown  in  Fig.  2;  the  distinctive  wineglass  shape 
had  not  developed  by  that  time  of  day. 


In  most  cases,  however,  the  meteorological  data  is  insufficiently 


DEPTH  (m) 


■omplete  or  precise  to  make  any  sort  of  theoretical  check  possible. 
\nother  difficulty  often  encounted  is  a  major  discrepancy  between  the 
jater  heat  rise  (profile  area)  and  the  cumulative  entering  heat  inferred 
:rom  the  meteorological  data;  it  is  suspected  that  shallow  temperature 
profiles  taken  close  to  the  observing  ship  can  be  substantially  distorted. 
■>ome  high  temperature  rises  under  especially  low  wind  conditions  have 
been  reported  by  Bruce  and  Firing  (1974)  (  }.6°C)  and  by  Romer  (  19691 
(4.4°C).  Although  these  data  sets  too  are  meteorologically  lacking,  it 
is  noteworthy  that  with  each  the  profiling  equipment  was  deployed  well 
away  from  the  main  experimental  platform. 
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1.2  Effect  of  transients  on  sonar  performance 

To  study  this, 

the  first  step  is  to  convert  the  transient  temperature  profile  over 
the  full  mixed-layer  depth  to  a  sound-velocity  profile.  Use  is  made  of 
a  sound-velocity  equation  such  as  that  of  Frye  and  Pugh  (1971).  From  the 
resulting  velocity  profile  may  be  calculated  the  transmission  loss  from 
source  (sonar  projector)  to  receiver  (target).  For  active  sonar  the 
loss  figure  in  dB  is  doubled  because  of  the  two-way  transmission  path. 
Service  documents  normally  discuss  the  topic  from  the  standpoint  of  ray 
theory.  More  precise  results  are  obtained  using  mode  theory,  as  explained 
typically  by  Pedersen  and  Gordon  (1965).  Mode  theory  takes  account  of  the 
sound  frequency  used,  and  ascribes  finite  transmission  losses  to  the 
"shadow  zones"  where  according  to  ray  theory  sound  does  not  penetrate  at 
all.  The  benefit  is  obtained  at  the  expense  of  extra  computing  complexity, 
however.  It  should  also  be  pointed  out  that  neither  theory  takes  account 
of  scatter  at  reflections  from  the  surface  due  to  roughness  of  the  latter. 

Fig.  3  shows  two  transients  with  corresponding  sound  velocity 
profiles.  These  profiles  are  idealized,  but  would  correspond  roughly 
to  conditions  of  moderate  solar  heating  and  windspeeds  of  5  and  10  knot. 

The  undisturbed  isothermal  profile  may  be  taken  to  represent  a  windspeed  of 
over  20  knot.  The  sound  velocities  at  the  source  depth  are  marked,  and  may 
be  compared  with  that  at  the  bottom  of  the  mixed  layer.  Some  hypothetical 
ranges  corresponding  to  these  profiles  are  shown  in  the  figure,  as  a 
function  of  receiver  depth.  Here  a  95dB  loss  from  source  to  receiver  is 
assumed,  and  mode  theory  is  used.  In  the  high-windspeed  or  isothermal  layer 
case  the  sound  velocity  throughout  the  layer  is  less  than  the  bottom  value 
so  that  the  duct  operates  to  full  advantage.  Taking  this  as  a  basis 
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for  comparison,  at  h)  knot  winJspevd  the  velocity  at  source  and 
neighbouring  levels  is  seen  to  he  raised  to  the  hotter,  value  approximate  iv. 
The  range  is  found  to  hi  somewhat  reduced,  more  particularly  for  the  shall¬ 
ower  receiver  depths.  At  3  knot  windspeed  the  source  velocity  is  seen  to 
be  weil  in  excess  of  the  bottom  value;  in  this  case  the  range  is  drasticall 
reduced,  from  20-30  km  to  about  7  km. 

Clearly  the  influence  of  transients  on  acoustic  transmission  is 
complex.  About  the  onLv  simple  generalization  to  emerge  is  that  a 
"break  point",  beyond  which  sonar  performance  is  seriously  impaired,  is 
reached  when  the  sound  velocity  at  the  source  depth  becomes  equal  to 
that  at  the  bottom  of  the  mixed  layer.  Fig.  4  shows  graphical Iv  the 
transient  temperature  rise  at  source  depth  4m  which  produces  this  equality, 
in  terms  of  the  depth  and  the  general  temperature  of  the  mixed  lavor; 
the  curves  are  calculated  using  the  sound  velocity  equation  referred  to 
earlier.  Hence  it  is  seen  that  the  surface  temperature  rises  cited  from 
Sverdrup  in  Section  1.!,  and  the  transients  shown  in  Fig.  1,  are  all 
capable  of  disrupting  sonar  performance  over  a  wide  range  of  mixed-layer 
conditions . 

The  transient  of  Fig.  2,  hv  contrast,  is  so  shallow  as  to  lie 
almost  wholly  above  the  source  depth.  Hence  it  should  serve  as  an  "over¬ 
head  mirror",  and  if  anything  improve  sonar  performance  by  eliminating 


surface  reflection  with  its  attendant  scatter. 
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1.3  The  present  studies 

Although  transient  observations  are  well  do.  umented,  very 
little  on  the  theoretical  side  is  to  he  found  in  the  literature. 
Accordingly  a  new  model  has  been  devised  to  enable  the  progressive 
development  ot  a  transient  during  the  day  to  be  calculated  frjm  the 
forcing  influences.  The  last  being  mainly  meteorological,  the 
prediction  of  sonar  impairment  hv  transient  then  becomes  an  extension 
of  meteorological  forecasting. 

The  model  is  successfully  fitted  to  a  limited  set  of  profiling 
and  meteorological  data  acquired  bv  the  author.  It  is  then  tested 
against  two  published  sets  of  empirical  rules  for  transient  formation, 
each  resting  upon  a  large  data  base,  and  there  is  found  to  be  broad 
agreement.  The  model  indicates  that  the  downward  transmission  of 
momentum  and  heat  due  to  turbulent  diffusion  is  shut  off  at  the  bottom 
of  the  transient,  completely  at  low  wind  speed  and  progressively  less  as 
the  speed  increases.  Hence  the  incoming  solar  heating  is  forced  to 
accumulate  above,  producing  the  "half-section  of  a  wineglass"  appearance 
characteristic  of  well-developed  transients.  Other  interesting  aspects 
of  transient  formation  arc  examined  also. 

The  model  is  used  to  predict  the  effects  of  some  important  practical 
environmental  changes  on  transient  formation.  The  results  of  a  typical 
afternoon  effect  prediction,  made  for  a  RAN  exercise  at  the  request  of  the 
Director  of  Oceanography  and  Meteorology,  are  given.  It  is  planned  that 
the  new  capability  be  added  to  the  AUSRAPS  sonar  range  prediction  scheme 
operated  by  the  RAN  Meteorological  Office  at  the  Naval  Air  Station,  Nowra 
NSW,  under  the  authority  of  the  Hvdrographer  RAN. 
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2.  The  model  used 

2.1  Diffusion  equations  for  heat  and  momentum 

These  are  based  on  the  well-known  empirical  concepts  of  eddy 
diffusivity  and  eddy  viscosity,  symbol  K,  dimensions  L2T-1;  has  not 
so  far  been  found  necessary  to  distinguish  between  them.  Salinity 
variation  has  not  yet  been  considered,  but  could  be  through  a  further 
equation  of  the  same  type.  The  diffusion  equations  are  frequently 
encountered  in  the  literature,  for  example  in  the  model  of  Mellor  and 
Durbin  (1975)  with  slightly  different  notation. 

The  heat  diffusion  equation  is  expressed  in  terms  of  temperature 

T(t,z) 

f  +  WO.  RZ(')  (1) 


where  t  is  time,  z  is  depth  measured  positively  from  the  surface,  and 
QR(t)  is  the  penetrative  heat  flux  (short-wave  radiation)  entering  the 
surface.  RZ(z)  is  the  penetration/depth  factor,  or  fraction  of  the 
entering  radiant  flux  absorbed  per  metre  at  depth  z;  it  is  also  the 
local  temperature  rise,  in  the  absence  of  heat  diffusion,  due  to  unit 
cumulative  flux. 

The  initial  condition  is  a  perfectly-mixed  layer  at  sunrise,  which 

provides  a  base  level  for  temperature  measurement:  T(0,z)  =  0  at  all  z. 

9T 

The  boundary  conditions  for  all  t  are:  z  =  0,  K  r— =QS(t);z=z, 

o  d  z  m 

3T 

—  =  0.  Here  QS(t)  is  the  total  non-penetrative  heat  flux  from  all  sources 

<3  Z 

entering  the  surface,  K  is  the  surface  value  of  K,  and  z  is  the  mixed-layer 

o  m 

depth , 


*  for  if  at  t  «  0,  W  -  W  at  all  z,  all  subsequent  W  are 
simply  increased  by  W  exp(-ift). 
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2.2  Diffusion  equation  for  turbulence 

The  most  widely  accepted  measure  of  turbulence  is  its  kinetic 
energy  E.  Making  the  usual  simplifications  its  local  rate  of  generation 
in  this  one-dimensional  model  is,  e.g.  Mellor  and  Durbin  (1975) 

ge  ‘  '  s6Q  -  e  (3> 


where  r  is  the  stress,  g  the  acceleration  due  to  gravity,  B  the  volume 
coefficient  of  expansion,  Q  the  heat  flux  diffused  downwards,  and  e  the 
dissipation  directly  to  heat.  If  turbulence  storage  and  transport  can 
be  neglected  then  G  =  0;  this  frequently-made  assumption  is  considered 
later.  However  since  the  present  situation  involves  substantial 
stratification  and  rapid  change  through  the  day,  we  wish  to  retain  these 
features  initially  at  least.  One  way  is  to  set  up  a  diffusion  equation 
for  turbulent  energy,  e.g.  Kundu  (1980,  equ.  5) 


3E 

3t 


_a 

3z 


(K 


3E 

3z 


)  +  G„ 


(4) 


The  disadvantage  here  is  that  the  additional  variable  E  is  introduced, 
so  that  a  final  closing  equation  is  still  required.  Some  modellers, 
e.g.  Kundu  (1980)  and  Marchuk  et  al  (1977)  add  yet  another  diffusion 
equation,  in  e  ,  making  a  further  closing  equation  necessary  still. 

The  author  instead  proposes  one  single  equation  for  closing 
(1)  and  (2),  based  entirely  upon  the  turbulent  diffusivity  K  (t,z). 

By  analogy  with  (4)  this  should  be  a  diffusion  equation  of  the  form 


3K 

3t 


_3 

3z 


«£> 


+  C„ 


(5) 


treating  K  as  a  scalar  property  which  can  be  stored,  diffused  and 
generated  -  in  fact  as  an  alternative  measure  of  turbulence  to  E. 
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These  are  the  principal  terms  in  a  model  developed  by  Nee  and  Kovasznav 
(1969)  for  such  conditions.  From  experimental  results  they  assign  the 
value  0.133  to  A.  This  value  turns  out  to  be  suitable  for  the  present 
work  and  is  adopted  henceforth. 


Initial  and  boundary  conditions  are  required  for  (9) .  The 
former,  at  sunrise,  follow  an  overnight  period  of  surface  cooling  and 
consequent  overturning  which  is  usually  sufficient  to  mix  the  upper 
layer  completely.  In  these  circumstances  K  is  assumed  to  commence  at 
some  value  K.  ,  uniform  throughout  the  mixed  layer  on  account  of  the 
overturning,  and  well  above  molecular  level  on  account  of  the  energy 
derived  from  both  the  overturning  and  the  overnight  wind  stirring. 

The  present  model  is  not  yet  sufficiently  advanced  to  determine 

from  overnight  conditions.  Instead  a  traditional  formula  is  adapted 
from  Sverdrup  (1942,  p.  494) 


K  (m2hr  !) 
w  ' 


0.41  W2 
0.05  W3 


whichever  is  less 


(10) 


where  W  (knot)  is  the  sunrise  wind  speed.  is  generally  made  equal 

to  K  ;  departure  from  this  is  explored  in  Section  4. 
w 

The  boundary  conditions  adopted  for  (9)  at  the  surface  and  at 
the  bottom  of  the  mixed  layer  are  both  —  =  0  ,  implying  zero  vertical 

flux  of  K  at  these  levels.  Here  the  equation  is  applied  by  expressing 
VO  and  BG  as  the  known  fluxes  of  momentum  and  heat  respectively, 
divided  by  K.  In  this  way  K  is  found  right  up  to  the  boundaries.  Improving 
the  precision  of  heat  and  momentum  flux  calculations  made  nearby. 

This  modelling  approach  is  less  radical  than  it  might  at  first 
seem.  If  turbulence  storage  and  transport  effects  were  small  enough  to 
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neglect,  then  (5)  would  reduce  to  G  =  0  just  as  (4)  would  to  G  =  0; 

k  L 

(6)  indicates  these  conditions  to  correspond.  We  are  in  effect  starting 

from  the  original  approximation  G  =  0  and  adding  empirical  terms 

E 

representing  turbulence  storage  and  transport;  the  unconventionality 

lies  in  basing  these  additions  upon  K  rather  than  E,  or  for  that  matter 

b  3  K 

upon  e.  As  terms  in  a  diffusion  equation,  —  and  7 —  (K— )  may  not 

dt  dZ  dZ 

3  E 

have  the  physical  standing  of  —  ,  but  they  do  not  seem  to  be  any  less 

3 1 

plausible  than  -r^-  (K— ) ,  ~  or  — ^ -  (K—S-) .  However  doubtful  the  new 
3z  3z  dt  3z  dz 

additions  may  be,  it  is  surely  better  to  take  some  account  of  turbulence 
storage  and  transport  rather  than  none  at  all. 
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2.3  Prescription  for  dissipation 

This  is  the  only  means  provided  for  fit! ing  the  model  to 
observed  data,  see  again  (9)  . 

Dissipation  thinking  for  the  present  model  was  influenced  hv 
early  experience  (Hill,  1980)  when  dissipation  was  omitted  altogether. 

A  reasonable  modelling  fit  was  obtained  at  low  wind  speed,  deteriorating 
at  higher  speeds  but  not  disast  rousl  y .  As  a  result,  present  praet  ic.  is 
to  make  the  dissipation  adjustable  from  zero  upwards.  Vhree  different 
prescriptions  are  examined,  each  with  its  multiplying  factor  of 
adjustment.  We  select  among  these  and  adjust  to  give  best  fit  of  the 
model  to  observat ional  data. 

Dissipation  is  commonly  considered  to  be  linked  to  turbulent 
kinetic  energy  througli  the  formula 

c  a  Ei/2/L  (10 

where  L  is  a  length  scale  provided  by  the  environment  .  i"  in  turn  is 
often  linked  to  K  through  another  formula 

K  <*  K  (12) 

Elimination  of  E  leads  to 

DIS  “  e/K  =  (K/lO)"  .  K  00 

where  we  introduce  a  numerical  "dissipation  factor"  V  to  be  adjusted 
later . 

As  an  alternative  to  (12),  Bradshaw  et  al  (1967)  replace  E 
in  (11)  by  the  stress  magnitude  |t|  .  This  leads  to 

DIS  <r  |  r  !  3'2/KL  =  (K/I.V-  .  VC3/-  .  F  I14) 

•  v 

A  third  possibility  is  simply  to  declare  the  turbulent  energy 
dissipation  rate  to  be  a  fixed  proportion  of  the  generation  rate. 
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Denoting  the  fixed  proportion  by  F 

DIS  =  (K.VC2.F  )/K  =  VC2.F  <15 ) 

However  if  turbulence  storage  and  diffusion  were  negligible  this  would 
equivalent  to  fixing  Ri,  which  is  theoretically  unrealistic.  Evidently 
these  three  prescriptions  for  DTS  are  members  of  the  one  family 

DIS  =  (K/L;)P  •  V<r~P  •  F  (16  1 

with  the  "dissipation  index"  P  taking  values  2,  ‘z  and  0  respect  ivelv; 
this  is  the  widest  range  possible  without  incurring  negative  indices. 

To  obtain  a  length-scale  L  we  follow  Mellor  and  Durbin  (19751, 
who  propose  the  ratio  of  the  first-  to  zeroth-order  moments  of  the 
F. 2  -  z  field.  Substituting  |t|  for  E  as  above  gives 


0  0 
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2.4  Method  of  solution 

Equations  (1),  (2),  and  (9)  are  solved  numerically  with  both  time 
and  depth  divided  into  discrete  intervals.  The  semi-implicit  (Crank  - 
Nicolson)  form  of  the  finite-difference  method  is  used,  as  explained 
typically  by  Razelos  (1973,  Section  4).  Each  of  (1),  (2)  and  (9) 
becomes  a  set  of  simultaneous  equations,  with  each  member  of  the  set 
representing  the  interaction  of  three  adjacent  deptn  elements  over  a 
single  time  interval.  The  set  is  solved  by  tridiagonal  Gaussian  elim¬ 
ination.  This  procedure  is  straightforward  for  the  linear  equations  (1) 
and  (2),  even  though  the  latter  involves  complex  algebra,  but  a  difficulty 
arises  with  (9)  on  account  of  its  non-linearity  in  K.  The  last  is 
overcome  by  expressing  each  K  (t,z)  as  the  sum  of  an  already-obtained 
approximate  solution  and  an  unknown,  supposedly  small,  correction;  the 
equations  are  reformulated  in  terms  of  the  latter  and  are  thereby  linear¬ 
ized.  The  resulting  corrections  are  applied  to  the  approximate  solutions, 
and  the  process  is  iterated  until  all  corrections  become  acceptably  small. 
Some  20  iterations  are  usually  found  necessary.  The  dissipation  scale  L 
is  calculated  in  the  same  iterative  loop. 

Immediately  following  the  solution  of  (1)  any  positive  temperature 
gradient  is  removed  by  "pooling"  the  least  possible  number  of  nearby 
elements.  This  is  to  simulate  the  natural  overturning  which  occurs  early 
and  late  in  the  day  when  surface  cooling  predominates.  Velocity  and 
diffusivity  are  pooled  also.  The  latter  is  increased  to  reflect  the 
turbulent  energy  gained  from  overturning,  using  an  adaptation  of  (6). 

At  present  the  standard  computing  timestep  is  5  min.  For  practical 
temperature  prediction  work  the  depth  elements  are  made  1  m  thick  to  20  m, 
thence  2  a  to  30  n,  5  m  to  50  m,  and  10  m  to  the  bottom  of  the  mixed  layer 
if  this  has  not  already  been  encountered.  However  for  systematic  studies 
the  upper  element  thicknesses  are  adapted  to  the  wind  speed.  At  the 
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lowest  speed  considered  here,  5  knot,  the  near-surface  detail  is 
found  to  require  \  m  elements  as  far  as  15  m.  At  the  highest,  15  knot, 
little  happens  close  to  the  surface  but  a  computational  instability 
often  encountered  in  this  type  of  numerical  analysis  is  found  to  enter; 
this  difficulty  is  avoided  by  amalgamating  all  of  the  1  n  elements  to  2  m. 

"Spot"  or  instantanous  data  values  are  entered,  usually  on  an 
hourly  basis.  The  program  interpolates  between  successive  values  to  suit 
each  intermediate  computing  timestep.  Two-hourly  outputs  are  usually 
sufficient.  The  solutions  are  found  to  be  well-behaved  and  not  to  require 
smoothing . 
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3.  Modelling  an  observed  transient 
3.1  Observational  data 

In  May  1976  the  research  ship  HMAS  l>  I  AMA.M  I  N.-\ ,  since  iln  unnissinncd, 
occupied  three  stations  in  the  Indian  Ocean  successive  1  ■'  tor  five  davs 
each.  Mixed-layer  temperature  profiles  were  taken  at  two-hourlv 
intervals,  cont inuous  ly  over  the  time  on  station,  using  standard  Sippiean 
expendable  bathythermographs  (XBT  s) .  All  necessary  meteoro logical 
observations  were  made  at  die  same  times. 

These  XBT  traces  suffer  from  a  slow  recorder  start,  causing  the  top 
5  m  of  the  profile  to  be  lost.  Moreover  the  warranted  absolute  accuracy 
is  0.2°C,  far  outside  that  needed  for  the  present  studies.  However  within 
any  one  record  the  level  of  coherence  is  considerably  better,  and  traces 
may  be  digitized  to  0.01°C  using  simple  magnifying  techniques  only.  To 
overcome  the  lack  of  an  absolute  temperature  baseline,  one  is  inferred  for 
each  profile  to  make  the  water  heating  (profile  area)  equal  to  the  net 
input  heating  from  all  causes,  cumulatively  from  sunrise  to  the  hour  .<f 
the  profile. 

fable  1  shows  the  general  conditions  obtaining  at  one  of  the  stations, 
at  location  12-13°S,  105-106°E.  Heating  was  only  moderate  on  account  of  the 
winter  season,  and  cooling  high  on  account  of  the  brisk  unceasing  winds. 

The  corresponding  temperature  profiles  are  shown  in  Fig.  5.  Those  for 
May  1  are  virtually  isothermal,  not  surprisingly  in  view  of  the  high  winds, 
and  are  accordingly  omitted.  On  May  2  there  was  heavy  overcast,  producing, 
temperature  rises  unsuitably  small  for  analysis.  The  May  3  profiles  are 
found  to  be  anomalous,  giving  peak  water  heats  much  higher  than  can  be 
reconciled  with  the  apparent  input  heating;  when  a  profile  is  moved 
horizontally  to  give  the  correct  included  area,  its  lowest  temperatures  a;  esc. 
to  go  substantially  negative.  This  unrealistic  behaviour  is  not  at  present 
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understood  .  However  on  May  -  and  3  .  end  i  t  ions  won  re  t avourab I e,  wltr. 
strong  heating;  these  profiles  lent1,  to  eoineide  lairiv  I  at  depth.  Si:i.v 
the  results  tor  the  two  Jays  are  similar,  t  iiev  ar<  •  •  •sited  it;  cr.i.r  : 
improve  aceuraev.  This  data  set  is  new  selected  ter  model  fitting. 

The  meteorologieal  observations  were  made  hv  trained  naval  personnel, 
and  comprised:  apparent  wind  speed  hv  hand-held  ane-vreter,  later  cor  retted 
to  true  speed  from  SATNAV  fixes  ot  ship's  position;  true  wind  direction 
inferred  from  wave  frtnr orientation;  dry-  and  wet-bulb  air  temperatures; 
sea  temperature  by  bucket;  a 'd  weather  state  generally,  including  cloud 
cover.  From  these  are  calculated  the  entering  radiant  heat  and  the  net 
long-wave  radiant  surface  cooling  using  charts  prepared,  by  lanes  (1966, 

Figs,  23  and  26),  also  the  net  air-sea  heat  exehangt  according  to  t he 
parameterisat ions  of  Friehe  and  Schmitt  (1976).  Fig.  6  shows  the 
penetrative  radiant  heating  and  the  overall  surfaec-coo! ing  rates,  also 
the  wind  speeds  and  directions,  all  composited  for  May  4  and  3.  Ihe  derived 
heat  unit  used,  C.m,  enters  naturally  witen  estimating  water  heat  Llie  area 
tinder  a  temperature  profile  scaled  in  °C  and  m:  I  C.n;  =  4.0  MJ  m  ". 

Data  is  also  required  on  water  clarity  for  calculating  the  depth 
distribution  of  heat  delivered  by  short-wave  radiation,  the  HZ  term  in  (1). 
Since  no  measurements  were  made  from  HMAS  DIAMANilNA,  the  optical  classification 
Oceanic  lb  is  estimated  from  a  regional  map  provided  hv  .lerlov  (1968,  Fig,  3(1). 
For  each  classification  Jerlov  tabulates  the  fraotioiv-oi  entering  radiant 
flux  which  remain  unabsorbed  at  particular  depths  (fable  XXI).  This  data 
for  classification  Oceanic  lb  is  shown  in  Table  2,  and  t  hat  for  the  middle- 
range  coastal  classification  Coastal  3  is  added  for  comparison.  For  the 
present  purpose  it  is  necessary  to  draw  a  smooth  curve  through  the  tabled 
data  points;  the  RZ  function  is  then  obtained  by  differentiating  with  respect 
to  depth.  Difficulty  was  experienced  at  first  with  manual  curve  drawing 
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[ABLE  2 

UNABSORBED  PERCENTAGE  OF  KADI  AX’!  FLL'X  AT  PARTICULAR  DEPTHS, 
FOR  REPRESENTATIVE  OCEAN T C  AND  COASTAL  TYPE  WATERS 

( AFTL-  'ERLOV  (1968)  TABLE  XXI) 


DEPTH 

(m) 

UNABSORBEl)  PERCENTAGE 

Oceanic  lb 

- r 

OF  RADIANT  FLUX 

Coastal  5 

0 

100 

100 

1 

42.9 

27.8 

2 

36.0 

16.4 

5 

25.8 

4.6 

10 

16.9 

0.69 

20 

0.02 

25 

7.7 

50 

1.8 

75 

0.42 

100 

i 

0.  10 

1 

i 

1 
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Radiant  penetration  -  depth  factors  for  representative 

oceanic  and  coastal  tvpe  waters. 
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leading  Co  roughnesses  in  the  RZ  curves,  leading  in  turn  to  anomalies 
in  the  graphed  model  results.  Hence  a  computerised  method  or  int erpo lat ion 
was  sought  instead.  The  standard  polvnominal  and  cubic  spline  methods 
can  produce  unrealistic  ("wiggly")  curves  when  the  datum  points  are 
sparse  as  in  this  case.  Accordingly  a  new  method  has  been  developed, 
giving  a  curve  more  acceptable  to  the  trained  eye,  still  with  continuity 
of  first  and  second  derivatives  throughout.  The  smooth  RZ  curves  then 
obtained  are  shown  in  Fig.  7.  This  development  is  of  some  general  interest 
and  is  proposed  for  external  publication  (Hill,  19R1). 


3.2  Fitting  the  modei  to  the  data 

The  foregoing  meteorological  data  lie  now  ente  ed  in  the  model  . 

For  a  first  trial  the  dissipation  inde  I’  is  set  at  • .  3  .uni  >'  at  a 
(no  dissipation),  4  and  lb  cessi ve i v .  lln  upper  >anol  of  Fig.  8 
shows  the  resulting  temperature  profiles  at  lines  1200,  1400,  1800  and 
1800.  These  profiles  are  farmed  hv  connecting  the  computed  points  with 
straight  lines  without  any  smoothing,  or  omission  of  points.  The 
processed  XBT  profiles  t  rom  Fig.  3  arc'  added  for  comparison. 

There  is  seen  to  be  broad  general  agreement  between  observation 
and  model,  improving  as  the  day  proceeds.  The  progressive  deepening  of 
the  transient  and  the  late-af tornoon  overturning  near  the  surface  are 
modelled  realist  leal  !v.  While  the  modelled  profiles  io  not  fit  i.v !  I 
enough  to  fix  F  precisely,  inspection  of  tl.eir  upper  regions  enables 
F  to  be  hounded;  clearly  zero  represents  too  little  dissipation  and 
lb  too  much.  We  provisionally  set  F  at  the  intermediate  level  4. 

The  other  two  dissipation  indices  give’  broadly  similar  results 
when  F  is  varied  through  an  appropriate  range.  The  lower  panel  of 
the  figure  shows  the  results  fur  all  three  indices,  each  with  F 
optimised  in  the  foregoing  manner.  Naturally  the  upper  parts  of  these 
curves  differ  only  slightly.  Beneath  the  transient  however,  in  what  we 
may  regard  as  the  "stem"  of  our  wineglass,  there  are  significant 
differences,  P  =  0  giving  the  most  heat  penetration  and  P  =  2  the  least. 
Indeed  with  the  latter  some  deep  negative  temperatures  built  up  in  the 
early  morning,  before  the  radiant  heating  becomes  significant,  remain 
undisturbed  all  day.  However  with  all  three  values  of  P  these  temperature 
rises  are  small.  Their  differences  with  the  "observed"  values  should  not 
be  taken  too  seriously  because  of  possible  errors  in  the  latter,  arising 
from  the  method  of  temperature  base  adjustment. 
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rhe  extreme  P-values  2  and  0  can  be  criticised  on  other  grounds. 

The  former  gives  noisy  results  near  the  surface  and  introduces  a  strong 
dependence  upon  the  doubtful  length  scale  L.  The  latter  gives  indifferent 
near-surface  performance  in  the  early  afternoon,  and  in  Section  4  is 
found  to  lead  to  computing  difficulties.  On  balance,  P  =  0.5  with  F  4 
seems  to  be  the  best  present  choice  among  the  three. 

Summarising,  we  have  a  range  of  models  to  choose  from,  all  broadly 
successful,  differing  only  in  the  prescription  adopted  for  turbulence 
dissipation.  The  implication  is  that  dissipation  plays  only  a  minor  part 
in  the  whole  process  of  transient  development.  We  shall  see  more  direct 
evidence  of  this  later,  and  meanwhile  pass  to  the  consideration  of  further 
input  data. 


4. 
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Testing  against  empirical  rules  for  transient  format  ion 

The  R.N.  Director  of  Naval  Oceanography  and  'Vteorolo;v  fDNOV, 
1968)  has  assembled  a  data  set  comprising  some  25000  bat  hvthermograph 
(BT)  soundings,  mostly  accompained  by  the  relevant  met  <  <>  --ul .  ig  i  nl 
observations.  The  BT  sWere  taken  during  1964-1967  from  l.N.  heather 
ships  on  Ocean  Weather  Stations  INDIA  and  UTIKi,  located  in  the  North- 
East  Atlantic  Ocean  at  latitudes  59°N  and  5Jh°N  respectively.  The  data 
was  analysed  by  DNOM  and  the  following  empirical  rules  for  transient 
behaviour  were  laid  down 

(1)  The  depth  of  formation,  to  the  top  of  the  transient,  is 
3,  6  or  9  m  for  a  mean  wind  speed  of  5,  10  or  1'  knot 
respectively,  over  the  three  hours  prior  to  furration. 

(2)  Signili'  nit  transients,  i.e.  those  clearly  visible  on 
BT  traces,  are  commonly  about  6  m  thick,  rarely  less 
than  3  m  and  not  often  more  than  9 

(3)  For  wind  speeds  over  1:5  knot  transient  format  ion  is 
unusual  . 

Another  data  set  has  been  assembled  hv  i.ibata,  Boston  and 
Boyce  (1965)  from  the  records  of  Ocean  Weather  Station  i\\PA  located 
in  the  North  Pacific  Ocean  at  latitude  50°N.  This  is  for  the  period 
1956-1959  inclusive,  and  for  the  summer  months  only  (April  through 
August)  when  most  t  ransients  form.  The  cited  data  compr:  ~s 

(i)  the  depth  to  the  top  of  the  shallowest  thermocline 
(i.e.  transient)  of  a  BT  always  taken  at  1620  local 
apparent  time,  and 

(ii)  the  average  of  5  wind  speeds  measured  over  the  preceding 


12  hours. 


TABLE  3 


TRANSIENT 

-DEPTH /WIND-SPEED 

RELATIONSHIP  AFTER 

TABATA  ET  AL  (1965 

,  FIG  3) 

WIND 

NO.  OF 

1 

!  DEPTH  TO  TOP 

OF  SHALLOWEST  THERMOCLINE  AT 

SPEED  , 

PLOT 

1620 

LOCAL  APPARENT  TIME 

(in) 

RANGE  | 

POINTS 

knot  | 

Range 

Mean 

Std.  deviation 

, 

4-6 

27 

0-10 

3 

2 

I 

(exc 1 .  2  point s) 

9-11 

1 

53 

0-17 

6 

2 

14-16 

38 

5-20 

12  j 

2 

I  1 

(excl .  1  point )  j  i 
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Any  depths  not  clearly  attributable  to  the  dav's  wind  speeds  are  rejected. 

The  results  are  reported  as  points  on  a  depth/wind-speed  linear  plot 
(Tabata,  Fig.  3)  extending  to  about  40  m  and  30  knot,  with  the  straight 
line  of  best  fit  superimposed.  The  last  is  not  used  here  because  the 
range  of  wind  speed  covered  is  wider  than  that  of  present  interest; 
instead  three  groups  of  plot  points  are  selected  in  the  vicinitv  of  the 
key  DNOM  wind  speeds,  viz.  5,  10  and  15  knot.  The  averace  depths  for 
the  three  groups  are  set  out  in  Table  3,  each  with  standard  deviation. 

They  are  seen  to  be  the  same  as  those  found  by  DNOM,  except  that  at 
15  knot  wind  speed  the  depth  is  13  m  instead  of  9  m. 

Thus  there  is  prospect  for  testing  the  model,  alre.idv  adjusted 
to  fit  data  for  a  narrow  wind  speed  range,  over  a  much  wider  range  even 
if  the  test  is  only  a  "broad-brush"  one.  To  this  end  a  sample  ocean 
location  is  taken  at  latitude  53°N  representing  all  the  Oh'S  -  cited. 

A  fine  midsummer  day  with  cloud  cover  r.ot  exceeding  two-eighths  is  assumed, 
and  the  remaining  model  inputs  are  based  on  climatological  data  for  the 
North-East  Atlantic  Ocean  region.  The  total  radiant  heating  for  the  day 
is  estimated  at  7.1  C.m*,  with  ncjn  peak  0.77  C.m.hr  ',  and  the  total 
surface  cooling  rate  at  (0.06  +  0.005  W)  C.m. hr  *  where  W  is  the  wind 
speed  in  knots.  The  general  mixed-layer  temperature  is  taken  as  15  C, 
corresponding  to  a  volume  coefficient  of  expansion  0.0002  C  As  before 
the  water  optical  classification  is  taken  to  he  Oceanic  lb.  Wind  speeds 
of  5,  10  and  15  knot  are  examined. 

This  package  of  idealised  meteorological  conditions  is  modelled 
as  before  with  P  *  0.5.  The  results  are  shown  in  Fig.  9,  in  a  separate 
row  for  each  wind  speed;  the  major  influence  of  the  last  upon  temperature 
rise  and  transient  depth  is  immediately  apparent.  Two-hourlv  temperature 
profiles  with  F  »  4  are  plotted  at  the  left-hand  side  of  each  row, with  temperature 


*  See  again  Section  3.1 
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and  depth  scales  chosen  to  give  all  profile  sets  a  broadly  similar  appearance. 
The  depth  to  the  top  of  the  transient  predicted  bv  the  DNOM  rules  is  add<  d  in 
broken  line.  To  make  a  test  the  intersection  of  the  DNOM  line  with  tile  1800 
profile  is  circled;  it  is  seen  to  he  reasonable-  close  to  the  near- isot  liermal 
"front"  of  the  profile  for  all  three  wind  speeds.  This  agreement  extends  to 
the  Tabata  ct  al  results  also. 

The  DNOM  rules  give  the  transient  thickness  simple  as  3-9  m;  the 
corresponding  range  of  transient  bottom  depths  is  added  to  each  diagram 
as  a  shaded  hand.  DNOM  temperature  reporting  is  to  0.1°c  so  that  a  rise 
of  this  amount  may  be  taken  as  defining  the  bottom  of  the  observed 
transient.  Accordingly  the  intersection  of  each  1800  profile  with  the  0. 1°C 
line  is  circled  also;  it  is  seen  that  these  intersections  lie  within  the 
bottom  depth  ranges  for  wind  speeds  of  3  and  10  knot,  although  not  for 
15  knot.  It  is  interesting  that  the  only  lack  of  consistency  should  be  at 
the  original  15-knot  level,  hut  here  wc  must  remember  that  the  whole  temporal  lire 
rise  is  only  twice  the  DNOM  reporting  threshold.  I  he  third  DNOM  rule,  viz. 
the  unlikelihood  of  transients  at  over  13  knot,  evidently  stems  from  their 
smallness  rather  than  their  absolute  non-occurrence. 

Thus  the  model,  originally  fitted  at  wind  speeds  of  10-13  knot,  is 
seen  to  respond  realistically  when  the  speed  is  reduced  to  5  knot.  We  can 
now  understand  why  sonar  operational  interest  is  1 imiled  to  the  wind  spec 
range  3-15  knot.  ,\t  highc-r  speeds  the  temperature  rises  arc  too  small  to 
matter;  at  lower  speeds  the  transient  is  so  shallow  as  to  lie  wholly  above  the 
sonar  projector,  thereby  providing  an  overhead  "mirror"  which  if  anything 
helps  sound  transmission  by  eliminating  surface  reflection  with  its  attendant 
scatter. 

The  central  column  of  temperature  profiles  in  the  figure  explores  the 
influence  of  F  .  At  10  and  15  knot  wind  speeds  it  is  much  as  in  Fig.  8, 
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the  speeds  being  similar.  However  at  5  knot  it  falls  to  a  level  of 
unimportance.  This  means  that  our  F  -value  original lv  chosen  for  high 
wind  speed  need  not  be  reconsidered  for  low  speed.  N i 1 -d  if f us  ion  prorilcS 
marked  'N'  are  added  to  the  diagrams  and  are  seen  to  coincide  with  the 
higher-F  profiles  beneath  the  transient. 

When  the  results  of  Fig.  9,  all  for  P  =  0.5.  are  repeated  with 
P  =  0  it  is  found  that  at  5-knot  wind  speed  the  iterative  solution  for  K 
(Section  2d)  fails  to  converge.  We  have  seen  that  the  sensitivitv  of  the 
modelled  temperature  profiles  to  F  decreases  with  decreasing  wind  speed, 
implying  that  the  dissipation  is  decreasing  also;  this  conclusion  is 
later  confirmed  in  Fig.  10.  Such  could  not  happen  with  the  P  =  0 
prescription,  which  ties  dissipation  to  VO  independently  of  wind  speed.  The 
implication  is  that  the  latter  nrescription  imposes  an  unrealistic  constraint 
upon  dissipation,  and  should  therefore  be  avoided. 

I'he  diffusivity  (K)  profiles  are  added  to  the  figure  for  later 
comparison  with  the  temperature  profiles,  but  discussion  of  them  is 
postponed  to  Section  5.2. 
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5.  Internal  behaviour  of  the  model 

5.1  The  gradient  parameters  VC-  ,  R<:  .md  P  1 S 

These  "bulk  system"  parameters  are  of  interest  ter  their  combined 
appearance  in  (9).  Since  gradients  are  traditionally  noisv  and  dinuult 
to  model,  we  do  this  using  the  idealized  data  of  Section  4  rather  than 
the  observed  data  of  Section  J.  I'hev  are  shown  in  Tig.  10  as  functions 
of  depth  for  the  same  conditions  as  Fig.  9,  all  in  the  unit  a”  .  Cluritv 
is  improved  by  plotting  the  snail  differences  (VC  -  HO  and  (VC  -  Ml'  -  M 
along  with  VG2  ,  with  ISO  and  DlS  appearing  as  curve  spacings. 

The  insensitivity  of  the  temperature  results  to  the  choice  of 
dissipation  prescription  has  already  been  mentioned,  so  the  comparatively 
low  levels  of  DIS  seen  in  the  figure  come  as  no  surprise.  ihe  even  smaller 
significance  of  dissipation  at  low  wind  speed  has  already  been  inferred 
from  the  central  profiles  of  Fig.  9. 

The  diagrams  show  VG2  and  BG  increasing  with  depth  to  a  peak  at  the 
bottom  of  the  transient.  This  peak  sharpens  as  the  day  progresses. 

Beneath  the  transient  they  collapse  abruptly,  VG'  to  zero  and  BG  to  a  low 
level  set  by  penetrative  radiation.  We  conclude  tha<-  the  downward  trans¬ 
mission  of  momentum  and  heat  due  to  turbulent  diffusion  is  largelv  supp¬ 
ressed  in  this  region,  for  all  three  wind  speeds  examined. 
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The  turbulent  diffusivitv  K 


Although  K  is  a  central  feature  of  the  present  work  we  must  not 
forget  that  it  is  an  empirical  concept  unsupported  by  present-dav 
turbulence  theory.  The  following  discussion  has  been  held  over  pending 
the  assembly  of  all  available  evidence. 

We  refer  again  to  the  modelled  results  of  Fig.  9.  Hie  range  ot 
K-values  encountered,  five  decades  in  all,  makes  the  use  of  a  logarithric 
scale  essential.  These  profiles  indicate  a  major  suppression  of  turbulent 
diffusivity,  developing  downwards  as  the  day  progresses.  At  any  part i i  tlar 
depth  within  the  transient  K  falls  from  early  morning  onwards,  reaches  i 
minimum  appropriate  to  that  depth,  and  changes  little  thereafter.  The 
upshot  at  1800  is  a  steeply-graded  K-profile  extending  from  the  surface, 
where  there  is  comparatively  little  change  from  the  assumed  sunrise  value, 
to  the  transient  bottom  where  K  has  fallen  through  several  decades.  Indeed 
at  5  knot  wind  speed  the  molecular  limit  is  reached  over  a  wide  depth  range, 
and  at  10  knot  it  is  closely  approached.  This  major  reduction  of  diffusivitv 
at  the  transient  bottom  evidently  acts  like  a  sheet  of  highly-insulating 
material  inserted  there.  The  downcoming  heat  cannot  penetrate  it  and  is 
forced  to  accumulate  above,  producing  the  characteristic  hal f-wi negl ass¬ 
shaped  temperature  profiles  shown  at  the  left  of  the  figure. 

It  is  possible,  in  principle  at  least,  to  infer  the  K  profiles 
directly  from  the  temperature  profiles  using  (1),  i.e.  without  modelling. 

Howe  and  Tait  (1969)  do  this  with  their  observed  temperature  profile  sets, 
of  which  a  single  profile  appears  in  Fig.  1.  They  are  able  to  follow  the 
fall  of  K  with  depth  as  far  as  a  factor  of  15,  beyond  which  the  calculation 
presumably  becomes  impracticable.  The  present  data  is  unf ort unate  1 v  not 
good  enough  for  any  such  analysis.  Nevertheless  it  seems  safe  to  conclude 
that  wineglass-shaped  temperature  profiles  such  as  those  seen  in  Figs.  1,2 
and  5  necessarily  connote  a  steep  fall  in  K  with  depth. 
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The  K-profiles  of  Kip.  9  all  exhibit  minima,  becoming  sharper  as  t  lie 
day  progresses,  which  delineate  the  transient  bottom  much  better  than  anv 
feature  of  the  temperature  profiles.  We  see  that  the  transient  deepens 
during  the  afternoon,  and  that  the  major  falls  in  K  occur  beneath  the 
transient  at  depths  where  the  latter  has  not  yet  arrived  Referring  again 
to  (9),  change  in  K  is  seen  to  be  controlled  by  the  group  (V-  -  BC  -  p;g) 
Of  these  terms  only  BC  persists  beneath  the  transient,  as  already  seen  in 
Fig.  10.  It  follows  that  the  falls  in  K  are  initially  caused  hv 
radiation  alone,  penetrating  to  sub-transient  levels. 

Within  the  transient  the  behaviour  of  K.  is  entirely  different.  The 
virtual  coincidence  of  each  set  of  profiles  in  this  region  shows  that  K 
remains  remarkably  constant  at  any  depth,  although  differing  widely 
between  depths.  The  group  (VC2  -  RC,  -  MS)  does  not  vanish  entirely  in 
the  region,  although  small,  ilowevi  r  (9)  shows  that  it  supplies  the  term 

4Z  (K-||)  as  well  as  |j:  ;  dividing  (9)  throughout  hv  K  makes  the 

relationship  clearer.  Fig.  11  shows  the  modified  terms  '  ■'  (K-  —  )  an.' 

K  '•  z  T  z 

1_  3K  .  -1 

K  ?lt  ’  1  the  unit  hr  ,  tor  two  wind  speeds  onlv  but  otherwise  for 

the  conditions  of  Fig.  H>:  we  see  that  within  the  transient  virtually  all 
of  the  drive  is  taken  up  by  the  diffusion  term,  leaving  the  storage  term 
negligible.  This  behaviour  is  consistent  wit  i  an  ad  hoc  assumption  made 
for  some  mixed-layer  models,  e.g.  that  of  Mel  lor  and  Durbin,  that  net 
generation  may  be  neglected  in  the  turbulent  energy  balance.  Such  an 
assumption  based  on  actual  ocean  observations  stands  by  itself,  but  here 
we  have  a  model  saying  the  same  thing  totally  unbidden.  Some  explanation 
is  called  for,  and  one  is  attempted  shortly. 

To  complete  the  discussion  of  K  we  recall  that  in  the  absence  of 
night-time  turbulence  analysis  it  is  necessary  to  assume  a  sunrise  value 
(lO  uniform  with  depth,  which  is  estimated  here  from  (10).  Fig.  12  shows 

the  effect  of  varving  by  factors 


of  0.25  and  4,  on  late  afternoon  temperature  profiles.  These  variations 
are  roughly  equivalent  to  halving  and  doubling,  respectively,  the  wind 
speed  during  the  foregoing  night.  We  see  that  the  effect  is  small  at 
15-knot  wind  speed,  ranging  down  to  negligible  at  5-knot.  Thus  bv 
late  afternoon  the  transient  evolves  to  a  state  largely  independent  of 
initial  conditions.  We  are  fortunate  that  ignorance  of  a  seemingly  kev 
input  parameter  does  not  prevent  us  from  modelling  transients  at  all 
wind  speeds  of  interest. 


WIND  STRESS  RADIANT  HEATING 

CORIOLIS  FORCE  SURFACE  COOLING 
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5.3  Feedback  loop  analogy 

We  now  attempt  the  promise  explanation  for  the  constancy  of  '< 
at  any  depth  within  the  transient 

!'ao  use  of  only  three  cqua  ions  in  the  model  makes  its  operation 
easier  to  understand.  Fig.  13  sh >ws  the  three  diffusion  equations  (I), 

(2)  and  (9)assembled  flowsheet -wi  ;e  in  a  closed  loop,  correspond inn  to 
the  actual  scheme  of  i omput at  ion .  We  start  with  (1)  and  (2)  at  the  top 
of  the  figure,  externally  forced  as  already  described.  K  plavs  the  same 
part  in  both,  so  that  its  influence  on  the  outputs  VC  and  BC  should  be 
similar.  However  these  enter  (9)  in  opposite  senses  and  in  differing 
degrees;  their  net  effect  after  processing  sets  the  rate  of  change  of  k. 

We  recognise  the  possibility  of  a  feedback  or  automatic-control  loop 
here,  tending  to  stabilise  K. 

To  check  polarities  of  changes,  consider  an  increase  in  wind  stress. 
This  increases  VG  as  it  enters  (9),  causing  a  progressive  increase  in  k. 
When  tills  is  fed  back  to  (2),  it  tends  to  "smooth  out"  the  velocitv 
profile  and  so  reduce  VG.  Hence  the  polarities  are  appropriate  for  stable 
control.  With  BC  however  the  reverse  applies  and  the  situation  is  unstabl 
since  BC  is  always  negative  the  tendency  must  always  be  for  K  to  fall. 

With  both  tendencies  in  operation,  stability  should  be  possible  as  long 
as  the  influence  of  VC  predominates.  Here  we  recall  that  VG  enters  (9) 
as  its  square,  but  BC  only  linearly.  Near  the  surface  where  VC  and  BC, 
are  high  the  possibility  of  control  is  most  apparent.  However  at  some 
depth  where  VG  and  BG  are  both  sufficiently  small,  the  relative  smallness 
of  VG2  should  cause  the  control  loop  to  collapse. 

Since  K  shows  every  sign  of  being  controlled  at  all  depths  within 
the  transient,  we  conclude  that  such  a  feedback  loop  is  operating;  or 
rather,  remembering  that  we  are  dealing  with  a  model,  that  transient 
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thermoc lines  behave  as  though  such  a  loop  were  operating.  1'his  analogue 
accounts  for  the  near-constancv  of  K  at  all  depths  within  the  transient. 
It  can  also  account  for  Lhe  general  stability  of  the  results  making  any 
form  of  smoothing  unnecessary,  and  for  their  insensitivity  to  the  assumed 
initial  value  of  K.  Beneath  the  transient,  however,  there  is  no  feedback 
loop  and  K  can  only  fall  as  already  noted. 


mrrfi^  pao s  eumc-not  fii>&> 
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6.  Behaviour  of  Richardson's  number 

Ri  BG/VG2  is  easily  obtained  from  the  results  of  Fig.  10,  and 
is  plotted  against  depth  in  Fig.  14  for  time  1600  and  all  other  conditions 
the  same.  Depth  is  normalised  by  expressing  it  in  ratio  to  transient  depth, 
defined  as  that  for  minimum  K.  Starting  from  a  small  negative  value  due 
to  surface  cooling,  Ri  becomes  increasingly  positive  with  depth,  steeply 
with  low  wind  speed  and  less  so  with  high,  approaching  1.0  slowly  near  the 
bottom  of  the  transient.  Thereafter  it  increases  rapidly,  since  VG2  falls 
to  zero  but  BG  does  not  on  account  of  penetrative  radiation. 

This  regime  differs  from  some  others  described  in  the  literature. 

Of  these  the  model  of  Meilor  and  Durbin  (1975),  henceforth  MD,  is  of 
especial  interest  because  of  some  close  parallels  with  the  present  work. 

It  uses  the  diffusion  equations  (1)  and  (2),  and  ignores  turbulent  energy 
storage  and  diffusion  so  that  (9)  reduces  to  (VG2-BG-DIS)  =  0.  The 
dissipation  prescription  (11)  is  used,  and  the  approach  to  calculating 
the  length  scale  L  is  that  adopted  here.  and  ,  the  differing  K's 

used  in  the  mechanical  and  heat  diffusion  equations,  are  calculated  from 
(12)  and  multiplied  in  each  case  by  an  algebraic  function  of  Ri.  These 
functions  are  such  that  both  and  fall  to  quasi-molecular  bottom 
limits  if  Ri  exceeds  0.25  approximately,  effectively  imposing  this  upper 
limit  on  Ri.  Such  a  limit  is  also  advocated  by  Turner  (1973)  from 
assembled  observations  of  incipient  turbulence. 

This  limitation  means,  in  the  present  notation,  that  BG  is  limited  to 
0.25  VG2.  From  the  MD  approximation  (VG2-BG-DIS)  «*  0  it  follows  that  DIS 
cannot  fall  below  0.75  VG2 .  Hence  the  turbulence  interplay  is  mainly 
between  generation  and  dissipation  directly  to  heat,  with  buoyancy  playing 
only  a  minor  part.  We  may  think  of  such  regimes  as  mildly  stratified,  in 
contrast  to  the  present  one  of  severe  stratification  where  buoyancy  plays 
a  major  part  and  dissipation  only  a  minor  one. 


To  explore  this  difference  further,  the  present  model  is  modified 
by  replacing  (9)  with  the  MI)  prescriptions,  and  then  tested  against  the 
same  data  sets  as  before.  The  resulting  1600  temperature  profile  for 
the  HMAS  DIAMANTINA  2-dav  composite  data  is  shown  in  Fig.  11,  with  the 
observed  and  modelled  (  F  =  4  )  profiles  repeated  from  Fig.  8  for 
comparison.  it  is  seen  that  the  MD  profile  does  not  lit  the  observed 
one  so  well,  being  too  shallow;  halving  or  doubling  the  MO  constant  of 
adjustment  a  from  its  recommended  value  0.1  makes  little  difference. 

The  Ri-profiles  are  shown  at  the  right,  confirming  that  Ri  readies  the 
values  0.23  and  1.0  near  the  bottom  of  the  transient  as  already  mentioned. 
The  results  for  the  idealized  OWS  data  are  shown  in  Fig.  16,  where  the 
outcome  is  similar.  To  put  the  comparison  beyond  doubt,  we  note  that  the 
two  Ri-curves  in  Fig.  15  run  roughly  in  the  ratio  0.23/1.0  over  much 
of  the  transient.  When  the  present  model  is  re-run  with  all  buoyancies 
increased  in  the  ratio  (0.23)  ',  tin.  MO  results  are  found  to  be  closely 
reproduced . 

On  the  basis  of  the  foregoing  evidence,  we  can  only  conclude  that 
the  natural  upper  limit  of  about  0.25,  which  is  sometimes  advocated  for 
Ri  in  the  mixed  layer,  does  not  apply  within  transients. 
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.  1*>.  I  emperat  ure  profiles  for  a  ir  l  «i  .urrvr  day  in  lat  iunir 
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7.  Practical  predictions  of  afternoon  effect 
7.1  Effects  of  key  environmental  •  hanges 

For  the  conditions  of  Fig.  9  (dissipation  index  !»  =  me  •  ■  .  , 

Fin.  17  shows  the  lf>OG  temperature  and  dilfusivitv  profiles  i\  r  res  peiu  i  in¬ 
to  the  following  input  variations. 

(SI  No  variation,  standard  case. 

(CC)  Cloud  cover  increased  from  two-eighths  maximum  to 

six-e igh thsapproximat ely .  Radiant  heating  halved, 
surface  cooling  unchanged. 

(TW)  Turbid  coastal  or  estuarine  water.  Coastal  t vpe  5 
instead  Oceanic  type  lb  (see  again  Fig.  ■'* )  . 

(WW)  Warm  tropical  water,  temperature  30°<;  instead  of 
1  5°C . 

(NO  No  Coriolis  acceleration,  location  equatorial 
instead  of  latitude  13°. 

The  effect  of  halving  the  radiant  heat  (profiles  CO  is 
naturally  to  reduce  the  surface  temperature  rise,  hut  the  reduction  is 
much  greater  than  simple  proportionality  would  predict.  This  is  because 
the  reduction  in  buoyancy  gradient  increases  turbulence,  and  so  enables 
the  diminished  heat  supply  to  penetrate  more  deeply.  At  critical  depths 
of  5  m  and  17  m  for  5-  and  10-  knot  wind  speed,  the  temperature  rises 
are  predicted  to  be  slightly  increased . 

The  influence  of  water  turbidity  (protiles  TW)  is  clearly  a  major 
one.  The  coastal  water  category  examined  here  is  the  central  one  of  five 
listed  by  Jerlov  (Table  XXI);  it  increases  the  surface  temperature  rise 
by  about  50%  at  either  wind  speed.  The  much-increased  radiant  absorption 
close  to  the  surface  effectively  removes  the  "stem"  from  t he  wineglass-shaped 
temperature  profile. 


$ 


Warmer  water  (profiles  UTi)  means  an  increase. 1  volume  coefficient 
of  thermal  expansion,  in  this  instance  from  .0001  to  .000 >  (  with  a 
proportional  increase  in  buoyancy  effects.  this  naturallv  makes  each 
temperature  profile  wider  and  shallower  leaving  the  radiant  stem 
unchanged . 

The  removal  of  Coriolis  acceleration  (profiles  NCI  favours  the 
transmission  of  momentum,  so  that,  conversely,  the  temperature  profiles 
become  narrower  and  deeper.  However  the  effect  is  comparatively  minor. 
Again  the  radiant  stem  is  unaffected;  the  temperature  profiles  S,  WW 
and  NC  are  all  seen  to  coincide  beneath  a  critical  depth  set  by  the  wind 
speed . 

Generally  it  can  be  said  that  the  profile  of  a  well-developed 
transient  represents  a  balance  between  the  two  main  forcing  agents,  sun 
and  wind.  Any  change  which  alters  this  balance  may  be  expected  to  have 
a  major  effect  on  the  profile. 
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7.2  A  typical  afternoon  effect  prediction 

The  following  gives  some  idea  of  what  can  l»e  offered  at  the  present 
stage  of  development  of  these  studies.  !n  1981  the  KAN  I'irecter  of 
Oceanography  and  Meteorology  requested  a  prediction  for  late  September, 
to  be  used  in  a  forthcoming  naval  exercise  near  Sydney  N.S.W.  (latitude 
3A°S)  .  The  basis  was  to  be  a  fine  sunny  day,  with  a  full  range  of 
windspeeds . 

Accordingly  the  representative  windspeeds  5,  7.  10  and  IS  knot  were 
selected,  and  the  radiant  heating  and  surface  cooling  during  the  dav 
calculated  as  for  Fig.  6.  The  water  optical  classification  was  taken  to 
be  Oceanic  lb .  The  modelled  temperature  profiles  for  local  time  1600  are 
shown  in  Fig.  18.  For  a  bulk  mixed-layer  temperature  of  17°C,  the 
mixed-layer  depths  giving  zero  sound-speed  excess  at  a  sonar  projector 
depth  of  A  m  are  calculated  from  Fig.  A  and  marked  upon  the  curves; 
these  are  the  minimum  Layer  depths  for  which  freedom  from  afternoon  effect 
can  be  assured. 
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7.3  Towards  a  routine  prediction  facility 

At  present  the  prediction  procedure  consists  of  three  separate 

parts 

(1)  A  computer  program  which  accepts  the  hourly  forecasts 
of  tiie  meteorological  variables,  and  finds  the  hourly 
radiant  heating  and  surface  cooling. 

(2)  The  main  program  which  accepts  the  above,  also  water 
optical  classification,  and  finds  the  hourly  mixed-layer 
temperature  profiles. 

(3)  Fig.  4  or  its  computational  equivalent,  converting  temperature 
profiles  to  sound-speed  profiles,  over  the  whole  mixed  layer. 
At  present  only  the  difference  between  speeds  at  projector 
depth  and  mixed- laver  depth  is  utilized. 

In  the  short  term  it  is  proposed  to  consolidate  (1)  -  (3)  above 
into  a  single  program.  A  compiled  version  will  be  filed  with  the 
CDC  7600  computer,  arranged  for  ready  access  by  the  RAN  Director  of 
Oceanography  and  Meteorology. 

In  the  longer  term  it  is  clearly  desirable  for  the  new  facility  to 
be  incorporated  in  the  existing  range  prediction  facility  AUSRAPS.  This 
question  has  not  yet  been  addressed. 
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8.  ('.inclusion 

What  do  we  learn  from  the  present  studies  of  afternoon  effort  ' 

Firstly  we  see  t  tat  it  is  basically  model  table.  Hu-  t  vpes  e< 
transient  profile  actually  observed,  with  their  , baraeterist ie  half- 
wineglass  shapes  and  strong  scaling  on  wind  speed,  are  rea 1 ist  i  a ' 1 v 
rendered.  This  is  achieved  using  only  standard  "tools  o;  the  t  r.idi  ", 
viz.  the  diffusion  equations  for  heat  and  momentum,  and  the  turbulent 
energy  balance,  here  ail  framed  in  terms  of  the  eddv-di f f us i vi t v  K. 

The  turbulent  energy  dissipation  prescript  iet.s used  are  also  adapted 
from  the  literature. 

Secondly,  on  using  dissipation  to  fit  t  lie  model  to  observation 
we  find  that  only  a  small  amount  is  needed.  Such  smallness  is  shown 
to  connote  severe  stratification,  with  gradient  Richardson's  number 
close  to  one,  over  most  of  the  transient  depth.  Since  this  finding 
contradicts  some  established  mixed-layer  thinking,  it  is  as  well  to 
review  the  evidence  briefly.  Of  the  HMAS  IHAMANT1NA  data,  that  for  two 
consecutive  sunny  days  is  composited,  thus  reducing  error  and  uncertainty. 
The  XBT  s  have  no  useable  temperature  base,  hut  when  individual  profiles 
of  the  diurnal  set  are  adjusted  for  heat  input ,  their  "stems"  are  found 
to  close  satisfactorily.  The  two  sets  of  empirical  rules  for  transient 
formation  are  "broadbrush",  but  are  mutually  consistent  and  rest  upon  large 
data  bases.  The  finding  of  small  dissipation,  made  using  either  data 
source,  is  almost  a  qualitative  one  and  therefore  not  easily  upset. 

By  contrast,  it  appeals  that  models  which  limit  buoyancy  formation  to  a 
low  level  cannot  fit  the  same  data  so  well. 

Some  interestinj  aspects  of  modelled  transient  behaviour  are  explored 
including  a  suggestion  of  feedback  loop  action  tending  to  stabilize  K.  Ms 
explored  are  the  effeits  of  varying  some  key  input  variables.  Generally 
it  may  be  said  that  t  ansients  represent  a  balance  between  the  contending 

influences  of  wind  and  sun.  Any  environmental  change  which  alters  this 
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balance  may  be  expected  to  have  a  major  effect  on  transient  development. 

On  the  utilization  side,  some  consolidation  of  computer  programs 
remains  to  be  done,  and  the  question  of  incorporating  the  whole  facility 
in  AlISRAPS  has  yet  to  be  addressed. 
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encountered  during  the  afternoons  and  evenings  of  calm  sunny  davs.  j 
It  is  caused  by  a  temperature  rise  near  the  top  of  the  otherwise  , 
isothermal  mixed  layer,  increasing  the  sound  speed  there  and  so 
refracting  the  sound  rays  downwards,  away  from  the  sonar  target. 

The  temperature  rise  profile  is  governed  by  wind  speed,  solar 
heating  and  other  meteorological  factors,  so  that  afternoon  effect 
prediction  should  become  an  extension  of  meteorological  forecasting. 
However  there  is  at  present  no  established  method  for  computing 
the  temperature  profile  from  the  meteorological  information. 

A  mixed-layer  model  is  presented  which  appears  to  meet  this  ! 
need,  subject  to  further  testing.  The  influences  of  some  kev 
environmental  factors  are  explored.  An  example  of  afternoon  effect  , 
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